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'~~P ! Abstract 

PLh ' We analyze a deterministic cellular automaton a = (ct" : n > 0) correspond- 

(-^ , ing to the zero-temperature case of Domany's stochastic Ising ferromagnet on the 

I hexagonal lattice H. The state space Sj^ = {—1,+!}^ consists of assignments of 

^ ' —1 or +1 to each site of H and the initial state cr" = {(T^ji^gn is chosen randomly 
I— with P{a^ = +1) = p € [0, 1]. The sites of H are partitioned in two sets A and B 

^ [ so that all the neighbors of a site x in A belong to B and vice versa, and the dis- 

^ ' Crete time dynamics is such that the a^'s with x ^ A (respectively, B) are updated 

'nI" . simultaneously at odd (resp., even) times, making a'^ agree with the majority of its 

, three neighbors. 

00 I In it was proved that there is a percolation transition at p = 1/2 in the 

' percolation models defined by a", for all times n G [l,oo]. In this paper, we study 

. the nature of that transition and prove that the critical exponents /3, and t] of the 

^ I dependent percolation models defined by o"", n € [1, oo], have the same values as for 

+^ ■ standard two-dimensional independent site percolation (on the triangular lattice). 

rl I Keywords: dependent percolation, critical exponents, universality, cellular automaton, 

■ zero-temperature dynamics. 
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1 Introduction 

The deterministic cellular automaton corresponding to the zero-temperature case of Do- 
many's stochastic Ising ferromagnet on the hexagonal lattice H j2j can be considered as 
a simplified version of a continuous time Markov process where an independent (rate 1) 
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Poisson clock is assigned to each site x G H, and the spin ax at site x is updated when 
the corresponding clock rings. The rule for updating the spin is to flip it if and only if it 
disagrees with two or three (a majority) of its neighbors. This model has been studied 
both rigorously and numerically in |,3, ; the numerical results about two critical exponents 
obtained there strongly suggest that the dependent percolation model defined by the 
limiting spin configuration cr°° is in the same universality class as ordinary independent 
percolation. 

In the zero-temperature Domany model studied here, the rule for updating the spins is 
unchanged, but the timing is different. This model has been previously used in numerical 
simulations and its dynamical as well as percolation properties have been studied 
in pp, [S] and [H]. In P], the existence of a percolation transition for the dependent 
percolation models associated with the spin configurations a" at time n, for all n G [1, cc], 
is proved; in 0, it is shown that for any n > 1, the crossing probabilities converge to 
Cardy's formula [7j when the lattice spacing 5 is sent to zero (the continuum scaling limit) , 
as in the case of independent percolation (at least on the triangular lattice T [^); in jH], 
the continuum scaling limit is analyzed in terms of cluster boundaries and is shown to be 
the same as for ordinary independent critical percolation on T. This last result strongly 
suggests that also the critical exponents defined at the critical point should be the same 
as for independent percolation. Indeed, due to the above mentioned result, one can use 
properties of the Stochastic Loewner Evolution [HI (which describes the scaling limit of 
ordinary critical percolation cluster interfaces) to compute certain critical exponents in the 
continuous model corresponding to the scaling limit. But unfortunately, the connection 
between critical exponents in the continuous and discrete models is not straightforward; 
more work and further results on the discrete model would be required to relate the 
discrete critical exponents with the continuous ones (see [TU]). 

Similar models on different lattices have been studied in various papers; see, for ex- 
ample, |llLll2H131ll4|ll5tlT61ll7j for models on 1/ and [TH] for a model on the homogeneous 
tree of degree three. Such models are also discussed extensively in the physics literature, 
usually on (see, for example, j2] and JH])- Their interest is tied to the fact that they 
can be obtained as zero-temperature limits of stochastic Ising models, a special class of 
Markov processes whose transition probabilities/rates are chosen so that the Gibbs mea- 
sures (for some Hamiltonian) at temperature T are invariant for the Markov process. In 
systems where there are multiple (infinite-volume) Gibbs measures for T below some crit- 
ical Tc, a subject of considerable interest is the t — oo behavior of the spin configuration 
cr* (with temperature T = Ti < T^) when the initial state is chosen from the (unique) 
Gibbs measure at T = T2 > T^. Studying the hmiting case where Ti = and T2 = 00 is 
the standard choice in much of the statistical physics literature (see, e.g., [201 )■ 

In this paper, we study the percolation properties of (say) -|-1 spins in the depen- 
dent percolation models, generated by the zero-temperature Domany cellular automaton, 
corresponding to ex"', with n G [l,cxo], when the initial state is chosen randomly with 
P((T° = +1) = p G [0, 1]. In [Ij, it was shown that there is a percolation transition with 
Pc = 1/2 for all values of G [1, 00]. We remark that is not 1/2 for n = 0, since the 
critical probability for independent site percolation on the hexagonal lattice is strictly 
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larger than 1/2, therefore the system is driven to criticahty by the dynamics, and this 
after just one time step. This is in contrast to what happens in the case of the continuous 
time model studied in [Hj, where it is believed that criticality is achieved again for p = 1/2, 
but only at time t = oo. Nonetheless, the nature of the percolation transition, in terms 
of critical exponents and/or scaling limits, is presumed to be the same in the different 
models. 

2 Definition of the Model and Preliminary Results 

In this section, we give a more detailed description of the model and present, for com- 
pleteness, results that were proved in jlj which motivate (and will be used in) the next 
section, where the main results of this paper are presented. 

Consider the homogeneous ferromagnet on the hexagonal lattice EI (embedded in 
so that the elementary cells are regular hexagons with side length 1 - see, for example. 
Figure with states denoted by o" = {crx}xem, = il? and with (formal) Hamiltonian 

{x<y) 

where denotes the sum over all pairs of neighbor sites, each pair counted once. We 

write A/bi(x) for the set of three neighbors of x, and indicate with 

cx^ay (2) 

the change in the Hamiltonian when the spin at site x is flipped (i.e., changes sign). 

Partition the sites of the hexagonal lattice IH into two subsets, A and B, in such a way 
that all three neighbors of a site x in ^ (resp., B) are in B (resp., .4). By joining two 
sites of A whenever they are next-nearest neighbors in the hexagonal lattice (two steps 
away from each other), we get a triangular lattice (the same with B - see Figure Q). The 
synchronous dynamics that we consider here is such that all the sites in A (resp., B) are 
updated simultaneously. 

We now define the cellular automaton cr", n G N, with state space Sm = { — 1, +1}^^, 
which is the zero temperature limit of a model of Domany ^.\, as follows: 

• The initial state is chosen from a Bernoulli product measure Pp, with -Pp(crQ = 
+ l)=p. 

• At odd times n = 1,3,..., the spins in A are updated according to the following 
rule: a^, x E A, is flipped if and only if A^.7i(cr) < 0. 

• At even times n = 2, 4, . . . , the spins in B are updated according to the same rule 
as for those in A. 
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We denote by cr°° the limiting state of the cellular automaton a^. a°° = limn^ooO""" 
exists with probability one J3] and, like cr" for 1 < n < oo, defines a dependent percolation 
model on H. 

The following observations are useful in understanding the behavior of the model and 
will help in the proof of our main result, Theorem El which is presented in the next section 
of the paper. 

• The values of the spins in A at time are irrelevant since after the first update 
those values are uniquely determined by the values of the spins in B. 

• Once the initial spin configuration in B is chosen, the dynamics is completely de- 
terministic. 

• A spin can no longer flip once it belongs to either a loop or a "barbell" of constant 
sign, where a barbell consists of two disjoint loops connected by a path. 

There is an alternative but equivalent way of describing the discrete time dynamics as 
a deterministic cellular automaton (with random initial state) on the triangular lattice T 
(corresponding to the set B of sites). Given some site x G T, group its six T- neighbors 

in three disjoint pairs {i/f,?/f}, {y^,yf}, {7/5, ?/g}, so that yf and are T-neighbors, 
and so on for the other two pairs. Translate this construction to all sites x G T, thus 
producing three pairs of sites {?/f,?/f}, {?/f,?/f}, {1/5,1/6} associated to each site x E T. 
(Note that this construction does not need to specify how T is embedded in M^.) Site 
X is updated at times m = 1,2,... according to the following rule: the spin at site x is 
changed from ax to —a^ if and only if at least two of its pairs of neighbors have all four 
sites with the same sign —ax- 

The models on the hexagonal and on the triangular lattice are related through a star- 
triangle transformation (see Figure [T] and, for example, p. 335 of [22] )• More precisely, 
the dynamics on the triangular lattice T is equivalent to the zero-temperature Domany 
dynamics on the hexagonal lattice EI when restricted to the sublattice B for even times 
n = 2m. 

To see this, start with T and construct an hexagonal lattice M.' by means of a star- 
triangle transformation such that a site is added at the center of each of the triangles 
[x, yf, ), {x, y^, 7/1), and {x, y^, y^). H' may be partitioned into two triangular sublattices 
A! and B' with B' = T. One can now see that the dynamics on T for m = 1,2,... and 
the zero-temperature Domany dynamics on H' restricted to B' for even times n = 2m are 
the same. 

An immediate consequence of this equivalence between the two cellular automata is 
that the two families of percolation models that they produce are also equivalent in an 
obvious way through a star-triangle transformation. To be more precise, the percolation 
models defined on T for times m = 1,2,... are the same as those defined on B by the 
zero-temperature Domany model for even times n = 2m. 

We now present (without proof) the results of jl] on the zero-temperature Domany 
model. Theorem^ says that the convergence to the limiting state is exponentially fast, 
while Theorem 121 identifies the critical point of the percolation transition, pc = 1/2. 
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Figure 1: A star-triangle transformation relating triangular and hexagonal lattice. 



Theorem 1. Let P^{n) denote the probability that a deterministic site in A flips after 
time n and similarly for Pein). Then, for any p, there is a constant c G (0, oo) such that 

P4H,Pf5(n) <e-^". (3) 

Theorem 2. Ifp > 1/2 (resp., < l/2j, there is percolation of +1 (resp., —1) spins in 
for any ri G [1, C)o] for almost every . 

Ifp =1/2: 

1. For n G [0, oo], there is no percolation in cr" of either +1 or —1 spins, for almost 
every . 

2. The mean cluster size is infinite in a" for all n E [1, C)o] (but finite for n = 0). 

3 Critical Exponents 

We will consider three percolation critical exponents, namely the exponents (3 (related 
to the percolation probability), z/ (related to the correlation length) and rj (related to 
the connectivity function). The existence of these exponents has been proved, and their 
predicted values confirmed rigorously, in recent papers PTIITU] , for the case of independent 
site percolation on the triangular lattice. Such exponents are believed to be universal for 
independent percolation in the sense that their value should depend only on the number 
of dimensions and not on the structure of the lattice or on the nature of the percolation 
model (e.g., whether it is site or bond percolation); that type of universality has not yet 
been proved. 



5 



Consider an independent percolation model with distribution Pp on a two-dimensional 
lattice L such that < Pc < 1- Let be the open cluster containing site x and IC^^I be its 
cardinality, then 6{p) = 6x{p) = Pp{\Cx\ = oo) is the percolation probability. Arguments 
from theoretical physics suggest that 9{p) behaves roughly like (p — Pc)^ as p approaches 
Pc from above. 

It is also believed that the connectivity function 

Tp{x,y) = Pp{x and y belong to the same cluster) (4) 

behaves, for the Euclidean length ||x — y|| large, like \\x — yW'"^ if p = pc, and like 
exp (— ||x — y\\/C,{p)) if < p < Pc, for some ^(p) satisfying ^(p) — oo as p t Pc- The 
correlation length ^{p) is defined by 

^{py^ = „ lim l-jr^ — n ^og Tp{x, y)\ . (5) 



\\x-y\\^oo I \\X - y\ 

C,{p) is expected to behave like {pc — p)^'^ as p '\ pc- 

It is not clear how strong one may expect such asymptotic relations to be (for more 
details about critical exponents and scaling theory in percolation, see [22] and references 
therein). For this reason, the logarithmic relation is usually employed. That means that 
the previous conjectures are usually stated in the following form: 

log^(p) 

lim 7—, r = f3, (6) 

plPa log[p - Pc) 

log Tp^{x,y) 

hm - — = -r], (7) 

\\x-y\\^^ log\\x - y\\ 

hm - — = -u. (8) 

pTpc log(pc - P) 

In the rest of this section and in the next one, 6x{p), Tp{x, y) and ^(p) will indicate the 
percolation probability, connectivity function and correlation length for independent site 
percolation on the triangular lattice T (identified with B). We will denote by 9x{p,n), 
Tp^n{x, y) and ^(p, n) the corresponding quantities for the percolation models on the hexag- 
onal lattice H at time n = 0, 1, 2, . . . . The main theorem of this paper is the following. 

Theorem 3. There exist constants < ci, 02,03,04 < 00 such that, \/n G [l,C)o], and 
X, y G EI and suitably chosen x', y' & B with | |x — x'| |, 1 1?/ — ?/'| | < 1, 

CiOx'ip) < ex{p,n) < C2 9xip), forpe (1/2, 1], (9) 
p'''Tpix',y')<Tp,nix,y)<p-''Tpix',y'), /or p G (0, 1/2], (10) 

e(p,n)=e(p), forpe (0,1/2]. (11) 

The next corollary is an immediate consequence of TheoremEland its main application; 
it says that the dependent percolation models defined by a", with n G [1, 00], are in the 
universality class of ordinary independent percolation. 
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Corollary 1. The critical exponents (3, rj and v exist for the dependent percolation mod- 
els on M defined by a"', with n G [l,c>o], and have the same numerical values as for 
independent site percolation on T. 

Remark 1. As already mentioned, the existence of the exponents (3, v andrj for indepen- 
dent site percolation on T has been recently proved, and their predicted values confirmed 
rigorously [F?lf7^ . 



4 Proofs 

Before we can prove the main results of this paper, we need some notation. We will 
denote by Pp^n the distribution of a" with initial density p of plus spins (i.e., with initial 
distribution Pp^ = Pp). 

For a site x G H, we will denote by M^ix) the set of its three neighbors in H. For a 
site a; G ;B, we will denote by Nb{x) the set of its six neighbors in B endowed with the 
graph structure of a triangular lattice (as explained in Section |2l- see also Figure [l}. 

We will call an W-path a path on the hexagonal lattice IH and a B-path a path on the 
triangular lattice B. A path whose spins are all plus (resp., minus) will be called a plus 
(resp., minus) path. Similarly, we will call an M-loop a (simple) loop on the hexagonal 
lattice EI and B-loop a (simple) loop on the triangular lattice B. A loop whose spins are 
all plus (resp., minus) will be called a plus (resp., minus) loop. Notice that constant- 
sign H-loops, doubly-infinite M-paths and "barbells" (two loops connected by a path) are 
stable for the dynamics, in the sense that, once formed, their spins will never flip again. 

IHI 

For two subsets C and D of EI, we indicate with {C < — > D} the event that some site 

IHI 

in C is connected to some site in Z) by a plus H-path, with {C < — > oo} the event that 
some site in C belongs to an infinite plus H-path. For two subsets C and D of B, endowed 

with the graph structure of a triangular lattice, we indicate with {C < — > D} the event 

that some site in C is connected to some site in D by a plus i3-path, with {C < — > oo} 
the event that some site in C belongs to an infinite plus i3-path. 

Proof of Theorem Let us first assume that x and y belong to B; we then take x' = 
X, = ?/ in Q - (fTUj) . The lower bound for 9x{p,n) in Eq. Q comes from the following 
bound 

6x{p,n) = Pp,„(x^^oo) (12) 

> Pp^o{{x < — > oo} n {x belongs to a plus H-loop}) (13) 

> Pp^o{x < — > oo) Pp,o(^ belongs to a plus H-loop) (14) 

> P^'Oxip) (15) 

> y^Up), (16) 

where we have used the fact that at time 1 (when the sites in A are updated for the first 
time) the dynamics transforms any constant sign i3-path into a constant sign H-path, the 
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FKG inequality and the fact that the events {0 < — > 00} and {0 belongs to a plus H-loop} 
are increasing (see, for example, |23); and the fact that the smallest H-loop contains 6 
sites. 

The lower bound for Tp „(x, ?/) in Eq. (jiup is obtained in a similar way. Mimicking the 
proof of the lower bound in Eq. 0, we have immediately 

rp,n{^, y) = Pp,n{x " y) > p^"^ Tp{x, y). (17) 

For the upper bound of Eq. Q we rely on the following observation. If no site in 
Mb{x) belongs to an infinite plus S-path at time 0, then, by the self- matching property of 
the triangular lattice, site x must be surrounded by a minus ;B-loop, which will produce 
a stable H-loop at time 1. Therefore, site x will not belong to an infinite plus H-path at 
any later time. Thus, 

O^p, n) < Pp,o{Mb{x) ^ 00). (18) 

Since the event {x < > 00} can be written as {cr° = +1} fl {A/b(x) < > 00}, using 

the FKG inequality we have 

Pp,o(x " 00) >pPp,oi^fBix) ^ 00). (19) 
From this and Eq. ()18|). we get 

eAP,n)<p-'e.,{p)<2 9,{p), (20) 

as required. 

For the upper bound of Eq. (fTUj). we first note that for bounded | |x — ?/| |, the inequahty 
is trivial by choosing C4 big enough so that the right-hand side of (llUj) exceeds 1. Next, for 

\\x — y\ \ large enough, we notice that unless {A/b(x) < — > N'siy)} at time 0, x and y must 
be separated by a minus i3-loop surrounding one of them or by a doubly-infinite minus 

IHI 

;B-path, and therefore it cannot be the case that {x < — > y} at any later time because at 
time 1 a stable minus H-loop or doubly-infinite H-path will be formed. This yields, for 
||x — ?/|| large enough, 

rp,n{^,y) < Ppfi{MB{x)^MB{y)). (21) 

Since the event {x < — ^ y} can be written as {cr° = cr° = +1} fl {A/b(x) < — > J\fs{y)}, 
using the FKG inequality we have 

Pp,o{x ^y)>p' PpfliMsix) ^ Msiy)). (22) 
From this and Eq. (pT|). we get 

rp,„(x, y) < p'"^ rp(x, y), (23) 

as required. 
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If X and y belong to A, the proof is analogous, but one has to consider shghtly different 
events. In this case we take x' G J^b{x) and then we have, for the lower bound in Eq. Q, 

6:^{p,n) = Pp,„(x oo) (24) 

> Ppfi{{x' < — > 00} n {x belongs to a plus H-loop}) (25) 

> Ppfl{x' < — > 00) Ppfl{x belongs to a plus H-loop) (26) 

> p'e^ip) (27) 

> Y^O^ip), (28) 

where ||a; — = 1 because x' is an H-neighbor of x. 
Analogously, for the lower bound in Eq. (fTUj) . we get 

Tp,n{x, y) = Pp,n{x <^ v) > p^^ Tp{x' , y') , (29) 

with x' G Mm{x) and y' G Mmiv)- 

For the upper bound of Eq. Q, we notice that, if no site in Uz£j^jj(x)-f^B{z) belongs to 
an infinite plus i3-path at time 0, then, by the self- matching property of the triangular 
lattice, site x must be surrounded by a minus S-loop that will produce a stable H-loop 
at time 1. Therefore, site x will not belong to an infinite plus H-path at any later time. 
Thus, 

0,ip,n) < Pp,o({U,eAre(.)A/B(2;) " oo). (30) 

Since for x' G Afm{x), {x' ^ oo} D {a° = +1, \fz G 7Ve(x)} H {U ,em^)M'B{z) " 
oo}, using the FKG inequality we have 

Pp,o{x' oo) > p^Pp,o(U2eA4i(x)A/'f3(z) « oo), (31) 

which yields 

0a.ip,n)<p-^e,ip)<8e,,{p), (32) 

with I |a; — x'l I = 1. 

The upper bound of Eq. (|TO|l . still trivial for bounded ||x — follows again from a 

similar observation for — ?/|| large enough: unless {l-iz€j\fM{x)J^Biz) ^ * ^z'eJ\fEi(y)-^Biz')} 

at time 0, x and y must be separated by a minus i3-loop surrounding one of them or by a 
doubly- infinite minus S-path. At time 1, a stable minus H-loop or doubly- infinite H-path 

IHI 

will be formed, making the event {x < — ^ y} impossible at any later time. This yields, for 
I |x — 1/1 1 large enough, 

rp,n{x,y) < Ppfl{VJz&Nn(x)^B{z) < > Uz'eMM{y)Wz'))- (33) 

Since for x' G Afmix) and y' G Mmiy), {x' y'} D {cr° = +1, Vz G J^mix)} (1 {cr°, = 

3 

+ 1, \fz' G Mmiy)} n {[Jzej\fM{x)J^Biz) < > Uz'eMiiy)-^Biz')}, using the FKG inequality we 

have 

Pp,oi^' < ' y') > Pp,o(U26A4i(x)A/'b(z) < > U2/6A4i(y)A/'B(/)), (34) 
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which yields 

Tp^nix, y) < Tp{x', y'), (35) 

with x' E A/h(x) and y' G Maiy)- 

The proof of Eq. (fTUI) in the remaining case (namely, x E A, y E B or vice versa) 
should now be clear. 

Eq. (fTTjl is an immediate consequence of Eq. (fTUj) and the definition of ^(p); it is 
enough to observe that 

M % 1 - || , HI [log ^pix', y') + C3 logp] I = (36) 

\\x-y\\^oo y \\x' - y'W J 

and 

M li^^ jT\\^ogTp{x\y')-Ci\ogp]\ =i{p)-^. □ (37) 

\\x-y\\^oo y \\x' - y'W J 

Proof of CorollaryU\ It follows from Eqs. (jHl) and pOj) that, for p G (1/2, 1] and ||x — ?/|| 
large enough, 

_ logci + log6'^/(p) ^ log 6'^(p,n) ^ _logC2jMogMp) /ggN 
log(p - 1/2) - log(p - 1/2) - log(p - 1/2) ' ^ ^ 

logri/2(x',?/') +C3log| logTi/2,n{x,y) \og Ti/2{x' , y') - C4 log | 

\ — n n — ~i — n rr — i — n n ■ y'^^) 

log\\x — y\\ log||x — yll log||x — yii 



Using these two equations, together with Eq. ()TTj) and the definitions of the critical 
exponents, and taking the appropriate limits gives the desired results. □ 



References 

[1] F. Camia, C. M. Newman, V. Sidoravicius, Approach to fixation for zero-temperature 
stochastic Ising models on the hexagonal lattice, in In and out of equilibrium: Proba- 
bility with a Physics Flavor (V. Sidoravicius ed.). Progress in Probability 51, 163-183, 
Birkhauser (2002). 

[2] E. Domany, Exact results for two- and three-dimensional Ising and Potts models, 
Phys. Rev. Lett. 52, 871-874 (1984). 

[3] C. D. Howard, C. M. Newman, The Percolation Transition for the Zero- Temperature 
Stochastic Ising Model on the Hexagonal Lattice, J. Stat. Phys. Ill, 57-72 (2003). 

[4] B. Nienhuis, private communication (2001). 

[5] F. Camia, C. M. Newman, V. Sidoravicius, Cardy's Formula for some Dependent 
Percolation Models, Bull. Brazilian Math. Soc. 33, 147-156 (2002). 



10 



[6] F. Camia, C. M. Newman, V. Sidoravicius, A Particular Bit of Universal- 
ity: Scaling Limits of Some Dependent Percolation Models, submitted, preprint 
|^rXiv:m atli.PR/03 08112| (2003). 

[7] J. L. Cardy, Critical percolation in finite geometries, J. Phys. A 25, L201-L206 
(1992). 

[8] S. Smirnov, Critical percolation in the plane: Conformal invariance, Cardy's formula, 
scaling limits, C. R. Acad. Sci. Pans 333, 239-244 (2001). 

[9] O. Schramm, Scaling limits of loop-erased random walks and uniform spanning trees, 
Israel J. Math. 118, 221-288 (2000). 

[10] S. Smirnov, W. Werner, Critical exponents for two-dimensional percolation. Math. 
Rev. Lett. 8, 279-744 (2001). 

[11] F. Camia, E. De Santis, C. M. Newman, Clusters and recurrence in the two- 
dimensional zero-temperature stochastic Ising model, Ann. Appl. Probab. 12, 565-580 
(2002). 

[12] L. R. Pontes, R. H. Schonmann, V. Sidoravicius, Stretched exponential fixation in 
stochastic Ising models at zero temperature, Commun. Math. Phys. 228, 495-518 
(2002). 

[13] A. Gandolfi, C. M. Newman, D. L. Stein, Zero-temperature dynamics of ±J spin 
glasses and related models, Commun. Math. Phys. 214, 373-387 (2000). 

[14] S. Nanda, C. M. Newman, D. L. Stein, Dynamics of Ising spin systems at zero tem- 
perature, in On Dobrushin's Way (from Probability Theory to Statistical Mechanics) 
(R. Minlos, S. Shlosman and Y. Suhov, eds.), AMS, Providence (2000). 

[15] C. M. Newman, D. L. Stein, Blocking and persistence in zero-temperature dynamics 
of homogeneous and disordered Ising models, Phys. Rev. Lett. 82, 3944-3947 (1999). 

[16] C. M. Newman, D. L. Stein, Equilibrium pure states and nonequilibrium chaos, J. 
Stat. Phys. 94, 709-722 (1999). 

[17] C. M. Newman, D. L. Stein, Zero-temperature dynamics of Ising spin systems fol- 
lowing a deep quench: results and open problems, Physica A 279, 156-168 (2000). 

[18] C. D. Howard, Zero-temperature Ising spin dynamics on the homogeneous tree of 
degree three, J. Appl. Probab. 37, 736-747 (2000). 

[19] J. L. Lebowitz, C. Maes, E. R. Speer, Statistical mechanics of probabilistic cellular 
automata, J. Stat. Phys. 59, 117-170 (1990). 

[20] A. J. Bray, Theory of phase ordering kinetics. Adv. Phys. 43, 565-580 (1994). 



11 



[21] G. F. Lawer, O. Schramm, W. Werner, One-arm exponent for critical 2D percolation, 
Electr. J. Probab. 7, Paper No. 2 (2002). 

[22] G. R. Grimmett, Percolation, second edition. Springer, Berlin (1999). 



12 



